Fourier Harmonics of High-p T Particles Probing the Fluctuating Inititial Condition 

Geometries in Heavy-Ion Collisions 
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Second Fourier harmonics of jet quenching have been thoroughly explored in the literature and 
shown to be sensitive to the underlying jet path-length dependence of energy loss and the differences 
between the mean eccentricity predicted by Glauber and CGC/KLN models of initial conditions. 
We compute the jet path-length dependence of energy-loss for higher azimuthal harmonics of jet- 
fragments in a generalized model of energy-loss for RHIC energies and find, however, that even 
the high-pT second moment is most sensitive to the poorly known early-time evolution during the 
first fm/c. Moreover, we demonstrate that higher-jet harmonics are remarkably insensitive to the 
initial conditions, while the different v n (N par t) vs. vi AA (N part ) correlations between the moments 
of monojet and dijet nuclear modifications factors remain a most sensitive probe to differentiate 
between Glauber and CGC/KLN initial state sQGP geometries. 

PACS numbers: 12.38.Mh,13.87.-a,24.85.+p,25.75.-q 



I. INTRODUCTION 

Heavy-ion collisions at the Relativistc Heavy Ion Col- 
lider (RHIC) indicate the production of an opaque (i.e. 
strongly jet-suppressing) pHl], strongly-coupled, fast- 
thermalizing medium that possibly needs to be described 
using methods derived from AdS/CFT string theory [6|. 
However, so far neither the initial conditions of the colli- 
sions nor the microscopic dynamics of the jet-energy loss 
are conclusively understood. 

Two models are commonly used to characterize the 
initial conditions. The Glauber model [7|, describing in- 
coherent superpositions of proton-proton collisions, and 
the "Color Glass Condensate" (CGC) @, given e.g. by 
the KLN model [9l-[ll| , where saturation effects are taken 
into account. They differ by their initial temperature gra- 
dients, their initial high-p^ parton distribution, and the 
distance travelled by each parton, leading to a different 
opacity estimate. In addition, both models exhibit large 
event- by- event fluctuations [ia - ila |. 

The jet-energy loss can either be described as multiple 
scatterings of the parton [i"6rl2lj |. specific for a weakly- 
coupled pQCD medium, or using the AdS/CFT corre- 
spondence where the problem of a parton stopped in a 
thermal medium is related to the problem of a string 
falling into a 5-dimensional black hole (22| - [25j . 

Experimentally, jet-energy loss is parametrized by the 
suppression factor Raa, defined as the ratio of jets pro- 
duced in A + A collisions to the expectation for jets being 
produced in p + p collisions 



Raa(pt) = 



dN A+A /dp T 
N coll dNP+P/d Pl 



(1) 



where N co u is the number of collisions, a theoretical pa- 
rameter (calculated within the Glauber model [l3j] ) de- 
pending on centrality (i.e., on the number of participants 

Np ar t). 



The first attempt to disentagle the initial condi- 
tions (Glauber vs. CGC) and the energy-loss mechanism 
(pQCD vs. AdS/CFT) while simultaneously describing 
the nuclear modification factor RAA(N part ) and the ellip- 
tic flow V2(N P art) for high-p^ particles was given in Refs. 
[2fj| - |28| . favoring CGC initial conditions and a strongly- 
coupled energy loss at RHIC. A similar ansatz was used 
in Ref. [29j j . although the Fourier components were not 
shown explicitly. 

To further study the differences of a pQCD and an 
AdS/CFT-like energy loss, we investigate the role of the 
energy dependence in the energy-loss prescription, ex- 
amine the power of the path-length dependence, and 
calculate higher-jet harmonics. We want to examine if 
a generic energy-loss ansatz that includes both a path- 
length and an energy dependence confirms the above 
conclusion that only CGC initial conditions and an 
AdS/CFT energy-loss can describe both the RAA(N part ) 
and the V2(N part ) appropriately. 

In the high-temperature limit, all dependences on the 
intrinsic scales of the system (T c ,Aqcd, etc.) disappear. 
Because of this, a generic energy- loss rate dE/dx is given 
by an arbitrary combination of dimensionful parameters 
constrained by the total dimension of the observable and 
the requirement that faster particles and hotter media 
result in a bigger suppression. We choose 



dE <- A \ 



-KP a T Z T Z - a+2 [x +h(cb)T,T}, (2) 



for an energy loss as a function of time r considering a 
jet starting at Xq and propagating in direction <fi with the 
coupling k. P is the momentum of the jet(s) considered, 
T is the temperature, and a, z are parameters controlling 
the jet energy (momentum) and path-length dependence, 
respectively. 

In the Bethe-Hcitlcr limit a = 1 and z = 0, while in 
the deep LPM pQCD limit a ~ and z ~ 1. If a = 
and z = 2, our model coincides with the model referred 
to as "AdS/CFT" in Refs. [H, HJ. However, on-shell 



2 



AdS/CFT calculations [2J-|25( show that a = 1/3 and 
z = 2 , thus we are going to consider a = 1/3 through- 
out the whole paper. However, one should keep in mind 
that a = 1/3 is a weak lower bound for falling strings 
used for illustration. We note that, as it is clearly shown 
in Refs. (22TI251]. z = 2 is only a lower limit correspond- 
ing to an "on-shell" quark whose stopping distance is 
If 1/T and whose initial energy is Eq ^> T. In a real- 
istic medium, the first assumption is likely to be violated, 
resulting in values possibly of z > 2. 

Please note that in contrast to Refs. (26l - [28j . n is a di- 
mensionless parameter. In a particular case of radiative- 
dominated scattering in the LPM regime, nT 3 ~ q [3Cj . 
In general, K measures the interaction cross-section, while 
T 3 is related to the density of scattering centers. This 
(like soft observables in general) can be used as a con- 
straint for k, but additionally one has to assume nearly 
complete thermalization on the timescale of jet propaga- 
tion throughout the system as well as a straight-forward 
relation between entropy density and multiplicity. While 
these assumptions are reasonable, they are not easily fal- 
sifiable, and, in particular, in the g-limit, one finds that 
k can not describe jet quenching and a realistic gluon 
density at the same time [30| . Therefore, we limit our- 
selves to fitting k to the most central data point of Raa 
measured at RHIC [27I . disregarding any interpretation 
in terms of the density of soft degrees of freedom, com- 
paring to the RHIC data on (xq, <ft)-avera ged Ra a versus 
centrality for a range of Ef ^6 — 9 GeV [271 ]. 

In a static medium, dE/dx ~ r 2 , while in a dynamic 
medium, dE/dx will aquire additional powers of r due to 
the dependence of temperature on t, implicitly included 
in Eq. ([2]). Here, we assume a ID Bjorken expansion [3l[ 



T(x,t) = T (x)(J±) 



1/3 



(3) 



with different values of r . Surprisingly enough, we find 
that even the high-p^ second Fourier moment (v% ) is most 
sensitive to the poorly known early-time evolution during 
the first fm/c. 

Considering that Raa is given by the ratio of jets in 
the QGP to jets in vacuum, one obtains the following 
formula for the nuclear modification factor from Eq. ([2]) 
for a jet starting at xq and propagating in direction (f> 



RAA(N par t,xa, 4>) = exp[-xO?o, 4>)] 



(4) 



with 



1 + a — n 
I- a 



In 



K 



1 ~ ■^l=£-f(£o><M,2) 



P, 



(5) 

where n is the spectral index (taken to be n ~ 6), K = 
k(1 — a), Pq is the jet's initial momentum and the line 
integral is 



I(x ,4>,a,z) = / t z T 



-a+2 



x + h(cf))T, r]di 



(6) 



In case of a = 1 (the Bethe-Heitler limit), 1/(1 — a) 
diverges, but since K = k(1 — a) — > 0, 



x(&o,<f>) = K(n-2)I(£ ,<t>,a=l,z). 



(J) 



From the above calculated RAA(N par t, xo,<p), the 
RAA{N P art) can be obtained by averaging over all possi- 
ble xq and (f> 



RAA(Np ar t) = 



d(j) J RAA{Npa r t,X ,4))TAA(xo)dxo 

2tt / T AA (x )dx 

^R AA (N pa rt,4>), (8) 

Z7T 



where the nuclear overlap function is used in case of the 
Glauber model and Taa = P 2 / Pq for the CGC prescrip- 
tion. 

After having fixed k, the v n (N part ) can be computed 
via 

f it v Jd(j} COS {n [0- lf; n }} R A a( N par t, (j)) 
V n {Np art ) = J - j— — — — . (9) 

J d(pRAA(Np ar t,(p) 

The Fourier density components e n and the reaction- 
plane axis ij) n are determined according to the initial den- 
sity distribution used in Ref. [HI 



and 



\ (r 2 cos(n(j))) 2 + (r 2 sin(n(/))) 2 
e n {t) = V _ (io) 

_ x (r 2 sin(n^))) 



i>n(t) = - tan" 

n (r z cos(n0)) 



(11) 



By definition, the impact parameter points into the x- 
direction, but event-by-event fluctuations introduce non- 
trivial and harmonic dependent ^> n 's. 

We checked that we reproduce the results of Ref. (2(| 
when our approach is simplified to the one used in this 
reference, where 

RAA{Np art ) = (e^ 7 ™), (12) 

and the line integral is 

dll™- 1 p(r + lv),m= 1,2,... (13) 



Here, m = 1 and m = 2 describe the path-length (I) de- 
pendence for pQCD and AdS/CFT- like energy-loss, rep- 
sectively. However, one of the main differences to our 
approach is that we also consider the energy dependence 
in the jet-energy loss. In other words, a = in Refs. [26l — 
[28j while we assume a = 1/3 [see Eq. ©] throughout the 
whole study. 

The above analysis can be extended to dijets, analo- 
gously to Ref. (28j . Dijet suppression is parametrized by 
the factor Iaa-, the ratio of the dijet suppression to the 
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FIG. 1: (Color online) V2 (top panels) and Raa (bottom pan- 
els) of high-momentum particles as a function of the number 
of participants, N par t, for z = 1 (left panel) and z = 2 [right 
panel, see Eq. @ for definition] and to = 1 fm. The RHIC 
data are taken from Ref. 12711. 



suppression of jets, defined experimentally as in [32fl and 
related to the parameters we defined earlier as in [28j | 



Iaa 



dN^/dpT 



- K (i*+i a ) 



RAAdN^JdpT ~ Raa 



(14) 



with the line integrals for the trigger (t) and away-side (a) 
jet. Please note that in contrast to Ref. (28[, our coupling 
k is equal for the trigger and the away-side jet. Clearly, 
a higher Iaa then Raa, as observed in Ref. [32|, implies 
that if one jet survives, the other one has a higher prob- 
ability to survive as well. In turn, this indicates that jets 
emitted from less dense periphery regions have a larger 
impact. Thus, as remarked in Ref. 32], comparing jet 
abundance to dijet abundance can be a sensitive medium 
probe. 



II. RESULTS AND DISCUSSION 
A. Mean Correlations 

We have carried out the procedure described in the 
previous section for a variety of impact parameters and 
initial conditions at RHIC energies. LHC energies have 
been mentioned in Ref. (33l . l3~ij and will be explored in 
detail in Ref. (35[. In this subsection, we focus on the 
mean correlations. The width of these correlations will 
be discussed in the following subsection. Please note that 
this width is the actual physical geometry fluctuation. 



FIG. 2: (Color online) V2 (top panels) and Raa (bottom pan- 
els) of high-momentum particles as in Fig. [T]for to = 0.01 fm. 



0.2 
0.15 

0.1 
0.05 





0.03 



KLN, z=2 — (b) 
KLN, z=5 - - {D > 
0.02 KLN, z=7 - - 



0.01 



-0.01 




0.1 0.2 0.3 0.4 0.5 

v 2 ' AA 



-0.03 -0.015 0.015 0.03 
v 4 Iaa 



FIG. 3: (Color online) Higher Fourier coefficients of dijet ob- 
servables, v\ AA vs. V2 (a) and v I 4 AA vs. va (b) for the Glauber 
model (reddish lines) and the KLN model (bluish lines) for 
z — 2, 5, 7 and to — 1 fm. A clear shift between the Glauber 
and the KLN model as well as a saturation effect for larger z 
can be seen. 



Throughout the paper, we distinguish four different 
cases. The Glauber model and the CGC/KLN pre- 
scription @-[l0|> both event- by- event and averaged over 
many events, i.e. in the latter case the initial conditions 
are smoothed over many events and the energy loss is 
calculated subsequently. 

In Fig. [T]we choose To = 1 fm, in line with recent hy- 
drodynamic calculations [36j . The figure shows that for 
both Glauber (red) and KLN (blue) initial conditions as 
well as for both pQCD-like [Fig. [2(a)] and AdS/CFT-like 
[Fig- LH C )] energy loss the RAA(N part ) can be described 
choosing an appropriate value for k. However, surpris- 
ingly enough, the results for V2(N part ) get close to the 
RHIC data when using KLN initial conditions for both 
pQCD-like [Fig. □]»] and AdS/CFT-like [Fig.Hfd)] en- 
ergy loss, while Glauber initial conditions underpredict 



4 



the data. Moreover, the difference between pQCD-like 
and AdS / CFT-like energy loss is rather weak. 

This result is a clear contradiction to the one shown 
in Ref. [13] and questions the conclusion that only 
CGC/KLN initial conditions and an AdS/CFT energy- 
loss can describe both the RAA{N part ) and the V2{N part ) 
appropriately. However, choosing a much smaller r , as 
done in Fig. reduces the absolute value of V2(N part ) 
for both pQCD and AdS/CFT- like energy loss scenarios, 
while it increases the difference between the pQCD and 
AdS/CFT results as seen in Refsjjll, [13] . Here, the dif- 
ference to the fit by Jia et al. [26| (black long dashed- 
dotted line), mainly seen for the pQCD-like energy- 
loss, is due to the additional energy loss dependence, 
parametrized by a = 1/3. Additionally, both plots show 
that there is a small discrepancy between the event-by- 
event and the averaged scenario 37']. 

We would like to mention here that in Refs. [26| dif- 
ferent values of tq were analyzed, nevertheless only the 
To = fm case was shown in Ref. [27J, leaving out a 
discussion about the physical meaning of To- 

Setting To = 1 fm means to assume that there is no 
energy loss within the first fm. PQCD does not give any 
excuse for this assumption and thus To = fm would 
be a natural assumption. However, To also describes 
the formation time of hydrodynamics which seems to be 
t ~ 1 fm 36]. On the other hand, setting To = 1 fm is 
also equivalent to the AdS/CFT result that the energy 
loss is suppressed at early times (due to the dE/dx ~ I 2 
dependence). Thus, it is important to note that the V2 
of high-px particles is sensitive to short-distance proper- 
ties, suggesting that there is either weak coupling with a 
To ~ 1 fm or strong coupling which in itself features the 
suppression of energy loss at early times. 

Extending the analysis to dijets and calculating the 
Iaa , our results meet the datapoint of Ref. [HI for both 
Glauber and KLN initial conditions when assuming that 
k is the same for both parts of the dijet, in contrast to 
Ref. 

Since the energy loss calculated by AdS/CFT [H-ii] 
clearly states that z = 2 is only a lower limit, this raises 
the question which power in the jet path-length depen- 
dece has to be included when an AdS/CFT-like energy 
loss is considered. A detailed discussion of this issue will 
follow in the next subsection (see the Fig. [5]). 

Before that, we will focus on correlations of higher 
Fourier coefficients of dijet observables, v\ AA vs. V2 and 
vl AA vs. U4. Fig. [3] displays those correlations for the 
Glauber (reddish lines) and the KLN model (bluish lines) 
for z = 2, 5, 7. 

A clear saturation effect occurs for larger z as well as 
a shift between the Glauber and the KLN model. Please 
note that while for the V2 correlation larger values of z 
of the Glauber model coincide with lower values of the 
KLN model, this is no longer true for the correlations. 
Therefore, the different means of the correlations between 
v n AA (pt) vs - Vu(pt) remain a most sensitive probe to dif- 
ferentiate between CGC/KLN and Glauber initial state 
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FIG. 4: (Color online) The fluctuation of V2 (top panels) and 
Raa (bottom panels) of high-momentum particles as a func- 
tion of the number of participants, N part , for z = 1 (left panel) 
and z = 2 (right panel) and to = 1 fm. The RHIC data are 
taken from Ref. [13]. 
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FIG. 5: (Color online) The fluctuation of the higher Fourier 
coefficients of dijet observables, v\ AA vs. vi (a) and v I i AA vs. 
va (b) for the Glauber model (reddish lines) and the KLN 
model (bluish lines) for z = 2, 7 and to = 1 fm. The blue 
area displays the region covered by the error bars of the KLN 
scenario for 2 = 2. The other cases show similar widths. 



geometries. 



B. The Shallowness of the Correlations 

After having focussed on the mean correlations in the 
last subsection, we now want to investigate the width of 
these correlations which is extremely important to con- 
clude about the significance of the previously shown dif- 
ferences in the path-length dependence of the energy loss 
and the initial states considered. 

Figs. H] and [5] are a repetition of Figs. [1] and G2 includ- 
ing the fluctuations for the event-by-event analysis. As 
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can be seen, those fluctuations are small for RAA(N par t), 
but the width for the v 2 (N part ) is rather large. Thus, 
in an experimental analysis it is less straighforward to 
distinguish between the different initial conditions than 
originally hoped for. Nevertheless, the CGC/KLN initial 
conditions seem to be favorable. 

In Fig. [5] it becomes obvious that the while larger 
Fourier coefficients become smaller, their width becomes 
larger, making it experimentally more difficult to disen- 
tangle the different scenarios. Please note that the scale 
of Fig. [SJb) is larger than the one of Fig. (3{b) , empha- 
sizing the very broad event- by- event fluctuations of the 
higher-moment correlations. Therefore, in order to draw 
an experimentally testable conclusion, it is necessary to 
always determine the mean and the width of the correla- 
tions considered. 

While the 2nd Fourier harmonics of jet quenching have 
been thoroughly investigated prl l28l . |39| , the sensitivity 
of higher harmonics has remained relatively unexplored. 
In Ref. [29j higher harmonics were mentioned but not 
shown explicitly. 

In a first step, we examine the differences in the ec- 
centricities between the Glauber and the KLN model, 
see Figure [HI Here it becomes clear again that while on 
average differences certainly exist (that are bigger for e-i 
than for 63), the width of the distributions is again rather 
large. This is true both for harmonics present on average 
as well as event-by-event, such as and harmonics only 
present once fluctuations are taken into account, such as 
e3, suggesting that higher harmonics of jet observables 
are rather insensitive to initial conditions. 

Fig. [7] displays the mean and the width of V3(N par t) 
and v^Npart) for Glauber and KLN initial conditions. 
While V3 is zero unless event- by- event fluctuations are 
taken into account, W4 is not too different between av- 
erage and event-by-event initial conditions. Please note 
that the absolute amount of our V3 seems to be a lit- 
tle larger than reported but not shown in Ref. In 
both cases, however, even a hypothetical ideal experi- 
ment capable of determining the impact parameter pre- 
cisely would be unable to distinguish between Glauber 
and CGC/KLN initial conditions using only jet harmon- 
ics, since event-by-event physical fluctuations are enough 
to drawn out model sensitivity. 

While higher-order coefficients are more sensitive to lo- 
cal gradients, they are also more susceptible to event-by- 
event fluctuations in initial conditions (hotspots, etc.), 
resulting in a larger 1^4 event- by- event fluctuation. 

Fig. [5] shows the sensitivity of Raa and v n to the mi- 
croscopic mechanism of energy loss, in particular to the 
power of the path-length dependence z. Here, we con- 
sider an impact parameter of b — 8 fm that maximizes 
v n . As can be seen, once Raa is fixed via the coefficient 
k, a residual sensitivity remains mostly in the Fourier 
components v 2 and V3. A saturation effect seems to oc- 
cur for larger values of z. 

Comparing the values of v 2 obtained by the Glauber 
and the KLN model in Fig. [5] to the RHIC data obtained 
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FIG. 6: (Color online) The eccentricity of the second (a) and 
third (b) Fourier components of the Glauber and KLN model, 
both event-by-event (red and blue lines) and averaged (red 
and blue dots). Since ez is only present when fluctuations are 
taken into account, it has to vanish for the averaged analysis. 



by the PHENIX experiment (black solid and dashed lines 



representing the datapoint at N, 



part 



1 25.7 and its error- 



bar) [27| show that KLN initial conditions seem to favor 
lower and Glauber intial conditions higher exponents of 
the path-length dependence. Note that z > 2 is allowed 
by AdS / CFT |22h25| , combined with a significant parton 
virtuality (which is reasonable since the typical stopping 
length is not If ^> 1/T). Since the virtuality can be di- 
rectly measured in jet-photon collisions, the dependence 
of the exponent on virtuality could become a quantita- 
tive signature of AdS/CFT dynamics. Thus, a simulta- 
neus measurement of v 2 , W3, and V4 could elucidate the 
microscopic mechanism of jet-energy loss. 

However, a question that naturally rises in this con- 
text is: Why are i> ra 's at high-p-r so insensitive to initial 
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FIG. 7: (Color online) V3 (top panels) and «4 (bottom panels) 
of high-momentum particles as a function of the number of 
participants, N pa rt, for z — 1 (left panel) and z = 2 (right 
panel) and To = 1 fm. 
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FIG. 8: (Color online) Raa (a) and v n (b-d) of high- 
momentum particles at an impact parameter of b = 8 fm as a 
function of the path- length exponent z, denned in Eq. ([2]) for 
ro = 1 fm. PQCD energy loss assumes that 2 = 1 [HL EaOill , 
for AdS/CFT on-shell partons z = 2, and AdS/CFT off-shell 
partons are z > 2 [22h25I ]. A saturation effect can be seen for 
all v n at large z. The RHIC data (black solid and dashed lines 
representing the datapoint at N par t = 125.7 and its errorbar) 
are taken from Ref. 12711. 



to reproduce the observed multiplicity distributions, and 
hence their (T) is similar, even if the local gradients of T 
might be very different. Therefore, if follows that hydro- 
dynamics and tomography lead to very different results. 

In conclusion, we investigated different Fourier har- 
monics of jet quenching at RHIC energies and showed 
that the second Fourier coefficients [v2(N part )] are re- 
markably sensitive to the initial time ro. If this tq = 
1 fm, as suggested by recent hydrodynamic calculations 
[36| . then the conclusion drawn in Ref. [27| that only 
CGC/KLN initial conditions and an AdS/CFT-like en- 
ergy loss can simultaneously describe the RAA(N par t) 
and V2(N P art) measured at RHIC can no longer be sus- 
tained. In contrast, for tq = 1 fm both pQCD-like and 
AdS/CFT-like energy loss seem to reproduce the RHIC 
data well if CGC/KLN initial conditions are taken into 
account. 

Moreover, we studied the microscopic mechanism of 
jet-energy loss by including an energy dependence and 
exploring the exponent z of the jet path-length depen- 
dence. We found that higher Fourier harmonics of jet 
quenching are remarkably insensitive to the differences 
between the Glauber and CGC/KLN model initial con- 
ditions. The different vf l AA (N part ) vs. v n (N part ) corre- 
lations between the moments of monojet and dijet nu- 
clear modifications factors remain a very sensitive probe 
to differentiate between Glauber and CGC/KLN initial 
conditions. 



conditions when hydrodynamics shows that the i? n 's of 
soft particles are extremely sensitive to initial conditions, 
leading to a systematic error of O (100%) in the viscosity 
[36j? In fact, the difference between these two regimes 
is not so surprising and can be readily understood physi- 
cally. While viscous forces are driven by local gradients in 
flow, jet absorption is driven by global differences in the 
integrated (-nP a T z T z ~ a+2 ) . The two effects are gener- 
ally not the same and can indeed be very different if the 
distributions (like initial distributions of energy density 
in a Lorentz-contracted nucleus) are not smooth. 

Both Glauber and KLN initial conditions are tuned 
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